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Abstract 


Crossflow instabilities dominate disturbance growth in the leading-edge 
region of swept wings. It is well known that streamwise vortices in a bound- 
ary layer strongly influence the behavior of other disturbances. Amplifica- 
tion of crossflow vortices near the leading edge produces a residual spanwise 
nonuniformity in the mid-chord regions where Tollmien-Schlichting (T-S) waves 
are strongly amplified. Should the T-S wave undergo double-exponential growth 
because of this effect, the usual transition prediction . methods would fail. 
Thus it is important to study interactions of this sort and to develop more 
realistic criteria for transition prediction. 

We modeled analytically the crossf low/Tollmien-Schlichting wave interac- 
tion as a secondary instability. The effects of compressibilty of the three- 
dimensional flow over a swept wing will be considered in a future analysis. 
The effects of suction are included, and different stability criteria are 
examined. The results are applied to Laminar Flow Control wings character- 
istic of energy-efficient aircraft designs. 

The work is an effort by personnel at Stanford University. At all times 
the work has been closely coordinated with an experimental program at Virginia 
Polytechnic Institute and State University (VPI & SU), which is examining the 
same problem. 
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1. Introduction 


1 .1 Motivation 

Recent research efforts in aircraft design are focused on the idea of 
Laminar Flow Control (LFC); see, for example, Pfenninger (1977b), Hefner and 
Bushnell (1977), and Bushnell and Tuttle (1979). The motivation for these 
efforts has been the realization that a turbulent boundary layer can account 
for up to 50% of the total drag on a vehicle during cruise. Maintaining a 
laminar boundary layer by delaying transition to higher Reynolds numbers and 
therefore increasing fuel efficiency is very promising and of practical im- 
portance to the industry. 

Transition to turbulence can be caused by many different things Morkovin, 
1969, 1977), but in the flight case transition occurs due to the unbounded 
growth of very small disturbances in the flow. It appears as though the most 
efficient means for maintaining laminar boundary layers is the use of very 
weak suction on the wing surface. The weak suction inhibits the growth of 
unstable disturbances in the boundary layer and delays transition. In the 
flight case this was demonstrated by Pfenninger and Groth (1961). The recent 
theoretical work of Reed and Nayf eh ( 1983) and experimental work of Reynolds 
and Saric (1983) and Saric and Reed (1983) demonstrate that, in the case of 
two-dimensional flows, the suction problem is well understood. The suction 
levels and distributions required to delay transition can be accurately pre- 
dicted . 

Modern LFC transport designs are expected to fly at subsonic speeds near 
the speed of sound. Their swept wings are being designed with advanced compu- 
ter codes (Bauer, Garabedian, and Korn, 1972; Kaups and Cebeci, 1977) and have 
supersonic flow on both the upper and lower surfaces. Typically, the upper 
surface of the two-dimensional airfoil is characterized by an extensive super- 
sonic flat-pressure region preceded by a leading-edge negative pressure peak 
and followed by a gradual shock-free recompression to subsonic flow with a 
subsequent rear pressure rise of the Stratford type (Pfenninger, 1977b); Alli- 
son and Dagenhart, 1978; Pfenninger, Reed and Dagenhart, 1980). 

Because of wing sweep, the flow over the wing is three-dimensional. The 
profile of the crossflow velocity component contains an inflection point. 
This situation is known to be dynamically unstable and creates streamwise 
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vortices that all rotate in the same sense Gray, 1972; Gregory, Stuart and 
Walker, 1955; Owen and Randall, 1952). These crossflow vortices are very 
important and dominate disturbance growth in the leading-edge region (Sro- 
kowski and Orszag, 1977). On the other hand, the mid-chord region is dominated 
by the classical Tollmien-Schlichting instability (T-S waves). Because the 
flow is locally supersonic, a complete three-dimensional T-S instability must 
be considered. 

One of the major unanswered questions regarding these instabilities con- 
cerns the interaction of crossflow vortices and T-S waves. Klebanoff et al. 
(1962) show that the onset of three-dimensionality is quickly followed by 
breakdown of the laminar flows. It is well known that streamwise vortices in 
a boundary layer strongly influence the behavior of other disturbances. Nay- 
feh (1981) shows that Gbrtler vortices produce a double-exponential growth of 
T-S waves. That is, the amplification rate increases exponentially in this 
case. Herbert and Morkovin (1980) show that the presence of T-S waves also 
produces a double-exponential growth of Gortler vortices, while Floryan and 
Saric (1980) show a similar behavior for streamwise vortices interacting with 
Gortler vortices. 

If crossflow vortices and Tollmien-Schlichting waves occur simultane- 
ously, it is possible that premature transition will occur (Pfenninger, 
1977a) . This possible interaction is especially important on the upper sur- 
face of supercritical LFC airfoils in the zero-pressure-gradient region. Vor- 
ticity developing in the leading-edge area may linger to the mid-chord region 
and affect stability there. Amplification of the crossflow vortices near the 
leading edge produces a residual spanwise nonuniformity in the mid-chord 
regions, where Tollmien-Schlichting waves are strongly amplified. Should the 
T-S wave undergo double-exponential growth because of this, the usual transi- 
tion prediction methods would fail. Thus it is important to study interac- 
tions of this sort and to develop more realistic criteria for transition 
prediction. 

1 .2 Review of Stability Theory 

Linear stability theory, or small disturbance theory, in which distur- 
bances are assumed in the form of linear, harmonic, traveling waves, plays an 
important role in the design of LFC systems. Using it, we can determine the 
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unstable, or amplifying, range of disturbance frequencies and amplification 
factors. The e^method (Smith and Gamberoni, 1956; Van Ingen, 1956; Jaffe, 
Okamura and Smith, 1970) and the modified e n and amplitude methods (Mack, 
1977) predict trends for changes in the mean flow that delay transition and 
are therefore very useful. However, they cannot be expected to predict the 
exact location of transition, because of the strong dependence of transition 
location on freestream turbulence levels. If the turbulence levels are low, 
the n-factor for predicting transition is about 15 (Runyan and George-Falvy , 
1979); for higher Reynolds numbers and turbulence levels, an n-factor of about 
ten has been suggested (Smith and Gamberoni, 1956; Van Ingen, 1956; Jaffe, 
Okamura and Smith, 1970). For LFC design, Hefner and Bushnell (1979) suggest 
a value of five. Saric and Nayfeh (1977) show how these values are changed 
when the nonparallelism of the basic flow is taken into account. 

For parallel, three-dimensional, incompressible flows, Gregory, Stuart 
and Walker (1955) derive the three-dimensional linear stability equations, 
including boundary-layer growth and streamline curvature. Then they determine 
a transformation reducing the three-dimensional temporal problem to a two- 
dimensional one. For flows over a rotating disk and a sweptback wing. Brown 
(1961) solves these equations numerically. Nayfeh and Padhye (1980) present a 
method for calculating neutral stability points for a flat-plate flow. Be- 
cause calculations of neutral stability (points separating stable and unstable 
flows) are extremely tedious and difficult in three dimensions, their work is 
significant. From their iterative scheme, they derive equations relating neu- 
tral and nonneutral disturbances. Cebeci and Stewartson (1980) identify an 
absolute neutral curve, called zarf, for the rotating disk. For given dimen- 
sionless frequency to, they find a neutral curve on which the growth rate is 
zero. That is, + p^ Can 4> = 0, y = x tan <p , where and p^ are the 
imaginary parts of the complex streamwise and spanwise wavenumbers, respec- 
tively, and $ is the direction of propagation. They then define the abso- 
lute neutral curve zarf on which both and p^ are zero. Padhye and 
Nayfeh (1981) do a nonparallel incompressible stability analysis on the X-21 
wing and find crossflow vortices to be the most amplified disturbance in 
regions of pressure gradient. In addition, they find the addition of nonpar- 
allelism to give more stable results in the front part of the wing and more 
unstable results in the rear pressure-rise region. 
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For parallel, three-dimensional, incompressible stability calculations, 
a computer code SALLY has been developed by Srokowski and Orszag (1977) that 
uses the e n -method for correlating the transition location. They calculate 
the maximum temporal amplification rate for a given dimensional frequency at 
each boundary-layer station from the parallel, incompressible stability equa- 
tions (the so-called envelope method). Then they use the real part of the 
group velocity to convert the temporal amplification rate into a spatial one 
(Gaster, 1962) and integrate along the path defined by the real part of the 
group velocity. In their limited calculations, they find stationary crossflow 
vortices to be the most amplified on LFC wings. Mack, in his spatial calcula- 
tions for the rotating disk (Mack, 1977) and Falkner-Skan-Cooke yawed wedges 
(Mack, 1978, 1979), also defines the direction of growth as that of the real 
part of the group velocity. In his work, he introduces his amplitude methods 
for correlating transition. Cebeci and Stewartson (1980) use the condition 
that da/dp be real, a condition also found by Nayfeh (1980a), to calculate 
an n-factor for the rotating disk. They start at a point of zarf and fix the 
group-velocity angle as a constant before marching. Nayfeh and Padhye (1979) 
establish a relation between three-dimensional temporal and spatial stabili- 
ties and a relation between spatial stabilities using the complex group veloc- 
ity. Malik (1980) and Malik and Orszag (1980) compare several methods of 
transition prediction using incompressible stability theory and conclude that 
the SALLY code is the most efficient. For the rotating disk, Malik, Wilkinson 
and Orszag (1981) then use SALLY to calculate temporal eigenvalues, which they 
convert to spatial eigenvalues using a group-velocity transformation. They 
then calculate ,n-f actors using the real part of the group velocity. Dagenhart 
(1981) considers stationary crossflow vortices and reports that his code MARIA 
adequately reproduces the stability results of the more complicated stability 
codes using the same physical disturbance model while using less than 2% of 
the computer time. 

Lekoudis (1979) and Mack (1979) evaluate the effects of compressibility 
on the stability of the boundary-layer flow over an infinite-span swept wing. 
They find that for T-S waves the inclusion of compressibility significantly 
reduces the maximum amplification rate and changes the most unstable wave's 
orientation. Mack finds for leading-edge-area crossflow that the incompres- 
sible theory overpredicts the amplification factor but has little influence on 
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the wave angle. El-Hady (1980) and Mack (1980, 1981) both report on the par- 
allel compressible stability of the flow over a 23°-swept wing with a super- 
critical airfoil shape. In their work, Malik and Orszag (1981) describe the 
computer code COSAL they have developed that efficiently computes temporal 
eigenvalues by finite differences. 

It is apparent that a correct three-dimensional stability analysis must 
include both compressibility and nonparallelism. Saric and Nayfeh (1975, 
1977) show that nonparallelism of the basic flow can be quite important under 
many circumstances. Nayfeh (1980a) uses the method of multiple scales (e.g., 
Nayfeh, 1980b) to formulate the three-dimensional problems but presents no 
numerical results. Nayfeh (1980c) then rigorously formulates the nonparallel 
3-D problem. He determines the partial differential equations governing vari- 
ations of the amplitude and complex wavenumbers and determines conditions on 
the group-velocity components, making the problem physically realistic. 
El-Hady (1980) presents some results on a 23°-swept wing with a supercritical 
airfoil shape. The work of Reed and Nayfeh (1982) represents the state of the 
art in three-dimensional calculations. Instead of jumping from wave to wave 
by locally calculating the most amplified disturbance in marching or specify- 
ing some artificial condition such as constant spanwise or chordwise wave- 
length, they formulate the problem by selecting a specific wave at some ini- 
tial point and then following that one wave along its trajectory, using the 
group-velocity ratio (Nayfeh, 1980a). A few results for the X-21 wing are 
given. They find again that the crossflow disturbances are the most amplified 
in the pressure-gradient region. 

If crossflow vortices and Tollmien-Schlichting waves occur simultane- 
ously, it is possible that premature transition will occur (Pfenninger, 
1977a). This possible interaction is especially important on the upper sur- 
face of supercritical LFC airfoils in the zero-pressure-gradient region. 
Vorticity developing in the leading-edge area may linger to the mid-chord 
region and affect stability there. Amplification of the crossflow vortices 
near the leading edge produces a residual spanwise nonuniformity in the mid- 
chord regions, where Tollmien-Schlichting waves are strongly amplified. Thus 
it is important to study interactions of this sort. 

For one to use the e n -method for the prediction of the transition loca- 
tion, it is necessary to know not only the respective growth rates but the 
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manner in which disturbances propagate in this three-dimensional interacting 
flow. Unfortunately, there is a dearth of analytical modeling and experi- 
mental data concerning the interaction problem. 

Two theoretical approaches are discussed in the next section to model the 
crossf low/T-S interaction (characteristic of swept wings) as a secondary in- 
stability. The first involves the weakly nonlinear combination of crossflow 
vortices with various oblique-wave solutions to the linear disturbance equa- 
tions. Craik's resonant triads (Craik, 1971) provide the foundation for this 
type of approach. The Craik model was considered for the onset of three- 
dimensionality of T-S waves in boundary layers. These types of models provide 
only limited understanding to interaction problems, because they specify one 
specific wavelength for secondary disturbances for a given frequency, provid- 
ing a very limited window. The secondary disturbances have been observed 
experimentally in a broad band of wavelengths (Klebanoff, 1962, Hama, 1960, 
Kachanov and Levchenko, 1982, Thomas and Saric, 1981, Saric and Thomas, 1983, 
Kozlov, Levchenko, and Saric, 1983). Moreover, interactions have been observed 
in flows where a Craik-type mechanism is invalid, e.g., subharmonic distur- 
bances in plane Poiseuille flow (Kozlov, 1982). 

The second approach described in the next section considers the T-S wave 
growth as a secondary instability in the presence of finite-amplitude cross- 
flow vortices. Herbert (1983) has successfully predicted the onset of three- 
dimensionality in plane Poiseuille flow and in flat-plate boundary layers 
using such an approach. Herbert considers the secondary instability of 
modulated finite-amplitude, two-dimensional T-S waves with respect to three- 
dimensional disturbances. For equilibrium plane Poiseuille flow, he solves a 
Hill-type system of equations with periodic coefficients and finds the differ- 
ent types of three-dimensional instability associated with different resonance 
conditions. His results (showing frequency and the broad spectrum of possible 
wavelengths as a function of amplitude and Reynolds number) compare very well 
with experimental observations. 

Boundary layers are more complicated, because they grow and do not exhi- 
bit the cross-stream symmetry of Poiseuille flow. Therefore, this analysis 
does not apply directly without some approximations. Herbert (1984) uses a 
shape assumption in which he neglects the nonlinear distortion of the two- 
dimensional T-S waves at finite amplitude and also ignores the weak variation 



of T-S wave amplitude on a viscous time scale. That is, he holds amplitude 
constant for his analysis so that his basic state (the flat-plate flow with 
two-dimensional T-S waves superposed) is periodic and a Hill-type system can 
then be solved for the interaction problem. The results of his parametrical 
excitation are consistent with available experimental data. 

1 .3 Objectives 

It is proposed to model analytically the crossf low/Tollmien-Schlichting 
wave interactions that are characteristic of flows over swept wings. This 
modeling is described in the next section. The effects of compressibility and 
nonparallelism of the basic flow will be included in a future analysis. The 
effects of suction in reducing the growth of these disturbances will be inves- 
tigated. Since the flow is complicated by the presence of crossflow vortices, 
the group-velocity-ratio criterion (Nayfeh, 1980a, 1980c) may not be the only 
applicable criterion. Thus other criteria will be eventually examined. 
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2 . Description of the Analysis 
2 . 1 Work in Progress — Analysis 

We consider a three-dimensional, incompressible stability analysis. For 
leading-edge-area crossflow. Mack (1979) finds that the incompressible theory 
overpredicts the amplification factor but has little influence on the wave 
angle. 

A. Basic State 

The spatial stability of the basic state (comprised of a three- 
dimensional boundary-layer flow over a convex or flat surface with a super- 
posed flow corresponding to steady streamwise vortices rotating in the same 
sense) is examined. This mean flow is typical of the upper surface of a 
transonic swept wing considered for Laminar Flow Control (LFC) in the flat- 
pressure region. This region is preceded by a leading-edge negative pressure 
peak where crossflow vortex instability is expected to be important. These 
vortices may linger to the zero-pressure-gradient region and affect stability. 

The three-dimensional boundary-layer flow is provided by the code of 
Kaups and Cebeci (1977), a code that solves the laminar compressible boundary- 
layer equations for a tapered, swept wing with wall mass transfer under the 
conical-flow assumptions. That is, the pressure gradient along the wing 
generators is assumed zero, and twist and tip and wing-body effects must be 
neglected. The governing equations are converted to a two-dimensional form by 
similarity transformations and then solved numerically by Keller's box method. 

Initially, we analyze the form of the solutions of the steady vortex 
flow. This was done by solutions of the disturbance equations for the vorti- 
ces. That is, we assume total flow quantities of the form 

A 

q(x,y, z,t) = Q g (x,y,z) + q(x,y,z) (1) 

Here, Q s (x,y,z) stands for the undisturbed three-dimensional flow provided 
by the Kaups-Cebeci code, and q(x,y,z) stands for the small, steady distur- 
bance due to the presence of the crossflow. Substituting the total flow quan- 
tities into the Navier-Stokes equations, subtracting the undisturbed part of 
this basic state (that is, the flow without the vortices), and linearizing, we 
find the first-order, nondimensional disturbance equations. We define the 
Reynolds number by 
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R = uV/v* (2) 

e e 

Here, U*, v*, and 6 * are the dimensional edge freestream chordwise 

velocity, kinematic viscosity, and displacement thickness, respectively. 

We confine our analysis to undisturbed flows that are only slightly non- 
parallel; that is, the normal velocity component V g i s small compared with 
the other components U s and W g • This in turn implies that all the undis- 
turbed variables must be weak functions of the streamwise and spanwise coor- 
dinates compared with the normal coordinate. In other words, 

U s = u s ^ x l»y» 2 l) » w s = W s (xi,y,zi) (3a) 

v s = e v J<*i. y. 2 i) ♦ v s = 0(i) (3b) 

We describe the relatively slow variations of the undisturbed quantities in 
the streamwise and spanwise directions by the slow scales xj - ex and zj = 
ez, respectively. Here e = 1/R. 

We describe the crossflow disturbance quantities in the form 

q(x,y,z;e, e y ) = e v lq 0 ( x i ,y,Zi> exp(18) +c.c.] (4) 

where c.c. stands for complex conjugate and 


50 

5x 


a v^ x l ’ Z i ) 


59 

5z 


^v^ x l ,z i 


(5) 


Here, e v describes the magnitude of the crossflow disturbance. The eigen- 
functions q Q are normalized so that their maximum value is one half. The 
quantities a y and p v are the complex dimensionless quasi-parallel compo- 
nents of the wavenumber in the chordwise and spanwise directions, respec- 
tively. Assuming the phase 0 to be continuously differentiable, we have 


5a v ^ dp v 
5z^ 5x^ 


( 6 ) 


For crossflow vortices, the wavenumber direction 4 , is deflected approxi- 
mately 90° from the potential-flow direction. 

If we adopt the convention 
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Z 01 ~ “ 0 » z 02 ~ ^0 » 


z 03 = v 0> z 04 = P0> z 05 = w 0> z 06 = DwO (?) 


we find that the equations governing the eigenfunctions of the crossflow dis- 
turbances can be written as a system of six linear first-order ordinary dif- 
ferential equations of the form 


6 


Dz 

on 



0 


m=l 




for n = 1, 2, 6. The a nm are the elements of a 6 

coefficient matrix whose nonzero elements are given in Appendix 
ary conditions for (8) become 


x 6 


A. 


( 8 ) 

variable- 
The bound- 


Z 01 = z 03 = z 05 = 0 at y = 0 (9a) 
Z 01 » z 03» z 05 * 0 as y > ® (9b) 

The system (8)-(9) constitutes an eigenvalue problem. Given the Reynolds 
number R and undisturbed mean-flow profiles, we determine the wavenumbers 
numerically. Since our problem is linear, it has six linearly independent 
solutions. The general solution is then a linear combination of these line- 
arly independent solutions, where the coefficients are determined from the 
boundary conditions. The procedure we use is to guess the eigenvalues and 
then numerically integrate the system (8) from the boundary conditions for y 
-* <*> to the wall at y = 0. In our analysis we employ the SUPORT code (Scott 
and Watts, 1977) and integrate with a Runge-Kutta-Fehlburg scheme using Gram- 
Schmidt orthonormalization to keep the solution vectors linearly independent. 
In general, our guess of the eigenvalues is wrong, and therefore one of the 
boundary conditions at y = 0 is left unsatisfied. We then use a Newton- 
Raphson scheme on this unsatisfied boundary condition to iterate and try to 
converge to the correct eigenvalue. The other conditions we satisfy are that 
the growth rate is a maximum and Wp/u) a is real. 

Once the correct eigenvalues are determined, we recover the eigensolu- 
tions and then construct the appropriate basic state for the interaction 
problem. 
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B. Crossf low/T-S Interactions Formulation 


In this section, two approaches are discussed to model the crossf low/T-S 
interaction as a secondary instability. The first approach models the T-S 
wave growth as a weakly nonlinear system and is similar to the work of Craik. 
The second approach considers the T-S wave growth as a secondary instability 
in the presence of finite-amplitude crossflow vortices. Herbert has success- 
fully predicted the onset of three-dimensionality in plane Poiseuille flow and 
in flat-plate boundary layers using this approach, known as parametric 
resonance . 

In both approaches, we superpose onto the basic state (that is, the 
three-dimensional boundary-layer flow over a convex or flat surface with a 
superposed flow corresponding to steady streamwise vortices, rotating in the 
same sense) two oblique, traveling, harmonic Tollmien-Schlichting waves. To 
this end, we assume small, unsteady perturbation quantities of the form 


e T q(x,y,z,t) (10) 

so that, considering (1), (3), and (4), the total flow quantities become 

As 

q(x,y,z,t) - Q s ( x 1 »7» z 1 ) + e v ( q o (x l » y ‘ z l> ex P (10) + c.c.) 

( 11 ) 

+ e T (q 1 (x,y,z,t) + q 2 (x,y,z,t)) 

The small dimensionless quantity e T represents the order of the amplitude of 
the two oblique Tollmien-Schlichting waves. Two oblique T-S waves are con- 
sidered to achieve a resonant or near-resonant condition. The Tollmien- 

Schlichting instability is assumed to be a secondary instability so that 

2 

is much smaller than e v and e. In other words, the terms of order are 

negligible in comparison with £j£y, and e^e . 

Then, using the method of multiple scales (e.g., Nayfeh, 1980b), we look 
for a uniform expansion of the q^-quantities (i = 1,2) in the form 

^(x.y.z.t) = f5 io ( x i >y« z i ,tj_) + e 9 11 (* 1 .y.Zj .tj) + ...‘ ] exp(10 1 ) (12) 

where 


ae 


i 


dx 


«iC x l .z^ 


(13a) 
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p i (x 1 , Zl ) 


(13b) 


de ± 

bz~ 


ae 


i 


at 


- U) 


(13c) 


The quantities a^, and oj are the streamwise wavenumber, spanwise 
wavenumber, and frequency of the two Tollmien-Schlichting waves, respectively. 
We substitute these quantities into the Navier-Stok.es equations and solve the 
resulting eigenvalue problem. 


Infinitesimal Crossflow Disturbance 

For growing boundary layers with infinitesimal crossflow disturbances, 
the parameter e v = 0(e) to account for both the effect of the vortices and 
the growth of the boundary layer at the same level. The disturbance equations 
describing the 0(e°) problem and the homogeneous terms of the O(e^) problem 
are identical in form to the crossflow disturbance equations given by (8) — (9) . 
(The basic state is again the three-dimensional boundary-layer flow over the 
wing without vortices.) Because of this, a solution for the inhomogeneous 
O(e^) problem exists only if a solvability condition is satisfied. The inho- 
mogeneous terms must be orthogonal to every solution of the adjoint homogen- 
eous problem. These solvability conditions depend on whether 

a, a a + a n and p » p + P~ (14) 

1 v 2 r l r v 2 

or not. If these resonance conditions are not satisfied, the solvability 
conditions yield two completely uncoupled equations describing the effect of 
non-parallelism on the two different amplitudes of the two T-S waves. If (14) 
is satisfied, we introduce detuning parameters and 02 defined by 

a v = a l “ °2 + eai 

(15) 

P, - h'h + t0 2 

and find the following coupled equations describing the amplitudes 

Aj (xf , zi , 1 1 ) and A 2 (xi,zj,tj) of the two T-S waves 
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dA 


g 


dA 


g. 


ii at s i2 ax 

d A 2 d A 2 


dA L 

+ S 13 dz L 

h A 

+ h 12 A 2 

exp(-ia 1 x 1 - 

(16a) 

dA 2 

+ g 23 dz l 

*2*2 

+ h 2l A l 

exp(-ia 1 x i - ic^z^) 

(16b) 


The resonance condition produces secular terms and causes the coupling of the 
amplitudes, leading to possible magnified growth of the T-S waves beyond that 
predicted by straight linear theory. The coefficients g^ t hj. , h^2» and 
1*2 1 are given in Appendix B. The quantities hi and h2 account for the 
nonparallelism of the mean flow, while h^ 2 an< i h 2 i account for the pres- 
ence of the crossflow vortices. 

Considering single-frequency, perfectly tuned disturbances, this system 
becomes : 

dA co fl dA. 

* + __i = h . A + h „ A„ 

dx a) . dz 11 12 2 

a 1 


co q a a 

— — + — ^ = h„, A, + 1 * 2^2 


dx co 0 dz 

a 2 


211 


(17) 


where (co„/<o and fco„/co are the group velocity ratios of the two T-S 

v P orl v p a'2 

waves. This system has the solution: 


A i " A i e 


^1 X+ ^2 Z 


A„ = A e 
2 2 


^1 X+ ^2 Z 


(18) 


locally, where the quantities A^ and A^ are the amplitudes of the two T-S 
waves as observed at the point (x,z) and 

<0 co 

p p 


1 


U 0 (0 

a 2 a 1 


a A 


co co 

Q A AAA O 

(h + h A /A ) — - fh + h A /A ) — 

v 1 12 2 l' co v 2 21 1 2 J ai 


‘I -» 


00 , 


Wn 


I "I 


CO 


a 1 


co 


a 2 


A A AAA AAA 

J h l ‘ h 2 ) + h l 2 Y A l ~ h 2l\ /A 2_ 


(19a) 


(19b) 


The quantities Si and S 2 are the growth rates associated with nonparallel- 
ism and the presence of the vortices. The total growth rates for the two 
waves are therefore 


= ~ Imag(a i ) - -J. Imag(P i ) + S l + 


P 


a i 


co L 

S 2 

a 1 


( 20 ) 
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Under these resonance conditions, we shall follow the T-S wave along its 
trajectory to ascertain its stability characteristics. Since the flow is 
complicated by the presence of streamwise vortices, marching along the charac- 
teristics defined by the group-velocity ratio may not be the only applicable 
criterion . 

Finite- Amplitude Crossflow Disturbance 

For finite-amplitude crossflow disturbances present in the growing 
boundary layer, we now assume that 0(e) < 0(e v ) < 0(1). Following the work 
of Herbert, the mean flow or basic state is described by: 

Q o (x,y,z,t) = Q g (x,y,z) + Eytq^Xj.y.z^ exp(i0) + cc] (21) 

where the eigenfunctions q Q have been normalized so that 

max | q | = 1/2 (22) 

0<y<» ° 

This expression describes a three-dimensional growing boundary-layer flow over 
a wing with embedded streamwise crossflow vortices. 

We neglect the nonlinear distortion of the crossflow vortices at finite 

amplitude, ignore the weak variation of crossflow wave amplitude on a viscous 

time scale, and apply Floquet theory. Under these conditions and considering 

the condition a, « a + a„ and B, => 8 + Bo to achieve resonance, the 
1 v 2 K 1 K v K 2 

first-order disturbance equations describing the two superposed oblique T-S 
waves are 


ia^u^ + Dv^ + 


= o 

(23 a) 

ia 2 U 2 + ° V 2 + 

1P 2 “ 2 

= 0 

(23b) 


/ a i + Pi \ 

+ ia l U s + ip l W s + — R— j U 1 + D Vl 

+ ia l P l ’ i “S + ie v | ( [ a 2 + °v] U 0 
+ P 2 "o) “2 + Pv U o“2j + e V V o Du 2 

+ e Du v = 0 (24a) 

V o 2 
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2 2 

, a 7 + \ 

(-iu + la 2 U s + ip 2 w s + R— -j “2 + DU s V 2 

+ia 2 P 2 " T ° 2|J 2 + i£ V|(_ a l " *vj “o 

+ Vo) U 1 ' Wi + Vo Du l 

+ e v I ^*o V l ' ° 

2 2 

/ + PA 

(- ico + iai U s + iP L W s + R ) V 1 + Dp l 

- A D 2 Vi + i £v (_^o + P 2 W ol V 2 
+ K u 2 + P v\] V o) + e V V o DV 2 

+ e v Dv o V 2 = 0 

2 2 

/ a 2 + P 2\ 

(- io> + i« 2 U s + 1P 2 « S + R j v 2 + ° P 2 

- j D 2 v 2 + l£ y (i_Vo + Vol V 1 
" _Vi + Vl] 7 o) + Vo D, l 

+ Ev CS? o v 1 - 0 

2 2 

/ a l + P l\ 

(- iu + ia L U s + ip x W g + g j + °Vl 

+ ^l p l ‘ 1 D S + i£ V (Vo“2 

+ |[ P 2 + P v] W o + Vo] “ 2 ) 

+ e V v 0 Dw 2 + e V DW o V 2 = ° 


(24 b) 


(25a) 


(25 b) 


(26a) 


L5 



(26b) 


2 L 2 

/ a 2 + h\ 

(- ^ + ia 2 U s + 1 P 2 W s + — — ) W 


+ ip 2 P 2 - j D 2 w 9 + ie 


2 + DW s V 2 


V (- W*! 
|[}l - P T ]» 0 + «!%{ »,) 


+ e„v Dw. + e,.Dw v = 0 

V o 1 V o 1 


The symbol indicates the complex conjugate. This coupled system of distur- 
bance equations is numerically solved by a spectral collocation method using 
Chebyshev polynomials. We apply an exponential mapping (Spalart, 1984) 


y = “ n o * n h (27) 

that transforms y = [0,®) to n = [1,0]. The value of n Q governs the 
concentration of points near the wall. Only odd Chebyshev polynomials are 
used in this formulation so that the infinite boundary conditions are automat- 
ically satisfied (Herbert, 1984). In this analysis, we use thirty polynomials 
giving a 123 x 123 real system of algebraic equations. The results are 

observed to be strongly dependent on the accuracy of resolution of the basic 
state. 

The spatial solution of this system provides the dispersion relation 

w = u>(a ,f3 ,R,£y) (28) 

The growth rates of the two T-S waves are therefore 

a, - - Imag(a.) Imag(p ) , i = 1, 2 (29) 

1 1 U) i 

a i 

Again we shall follow the T-S wave along its trajectory to ascertain its 
stability characteristics. The amplification factor n is then the integral 
along the trajectory of the growth rate. 


2 .2 Work in Progress-Calculations 

At present the interaction codes are being developed using X-21 wing 
data. This wing, derived from the NACA 65A210 airfoil, has leading and 
trailing edge sweepback angles of 33.2° and 19.1°, respectively. The x,y,z 
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coordinate system used in the stability analysis is shown in Fig. 1. The 
input upper surface pressure, suction, and calculated displacement thickness 
Reynolds number distributions are shown plotted against percent chord in Figs. 
2, 3, and 4, respectively. These quantities apply for a freestream speed of 
774.4 feet per second and a 14.66-foot chord. 


A. Crossflow 


The crossflow instability is expected to dominate in the pressure- 
gradient regions. The crossf low/T-S interaction we examine will occur at the 
leading-edge negative pressure-gradient region. Amplification of crossflow 
disturbances there produces a residual spanwise nonuniformity that could 
linger to regions where T-S waves are strongly amplified. For this reason, we 
focus our attention on the leading-edge area. 

Figure 5 shows typical streamwise and crossflow velocity profile distri- 
butions for the leading-edge area; these are computed at 1.5% chord. In this 
plot, 

* 

U s = dimensional chordwise speed 

* 

U = dimensional chordwise edge speed 
e 
* 

W s = dimensional spanwise speed 
* 

W g = dimensional spanwise edge speed 

Note the inflection point in the crossflow profile, a situation known to be 
dynamically unstable. Also of interest is the relatively small amplitude of 
the crossflow velocity in comparison with the streamwise velocity. The 
distribution of maximum (in magnitude) dimensional crossflow speed W max 
versus percent chord is shown in Fig. 6; the values are overall very small. 


Two dimensionless parameters commonly used both in experiment and theory 
to correlate crossflow-vortex appearance with freestream and boundary-layer 
conditions are the local crossflow Reynolds number 


and local shape factor 


Re 


CV 


W v 
max^lO 
* 

v 

e 



y 

J max 
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Distributions of these plotted versus percent chord are shown in Figs. 7 and 
8, respectively. Vortices appear when crossflow Reynolds numbers are on the 
order of 100 and shape factors are on the order of 0.3; these figures are 
consistent with our findings. Here y max t * ie l° ca l normal-direction 
location of ^ max and y^Q t * ie norma l - ditection location above y max 
where the crossflow speed is 10% of W max « 

A spatial stability analysis on this flow yields the following results : 
Figure 9 shows the amplitude distributions of two perpendicular disturbance 
profiles, one in the direciton of the wave angle, the other normal to this 
direction. These are computed at 1.5% chord. The larger amplitude portion 
lies in the direction perpendicular to the wave and is normalized to 1/2 to 
give proper meaning when an amplitude ey Is introduced. 

Figure 10 shows parallel and nonparallel calculations for the most 
unstable disturbance at the leading edge. For Padhye and Nayfeh (1981) 
nonparallelism is a stabilizing influence. Here we see that nonparallelism is 
initially stabilizing, then destabilizing. Figures 11, 12 show how that 
disturbance's wave angle and group velocity vary with percent chord. It is 
apparent that the most unstable disturbance lies in the direction perpendic- 
ular to the freestream and is a crossflow-type growth. The group-velocity 
direction is almost the same as the potential-flow direction. 

B . Crossf low/T-S Interaction — Infinitesimal Disturbance 

An appropriate basic state is formed and the interaction problems now 
considered. Because the Craik-type approach is expected to provide only 
limited insight into the interaction problem, only a few results are presented 
here to demonstrate the model. The model depends on knowledge of the initial 
conditions (that is, the amplitude ratio) of the two interacting T-S waves, as 
well as the amplitude of the crossflow disturbance. Moreover, the model 
provides a single disturbance wavelength for a given frequency. 

Figure 13 shows growth rate (defined by Eq. (20)) plotted versus T-S wave 
amplitude ratio for three different crossflow velocity magnitudes (0%, 0.6%, 
1.2%) for a chord location of 0.5% and a frequency of 325 Hz. As is evident 
from the plot, the stability of the oblique T-S wave is strongly dependent on 
the amplitude ratio. Predicting this ratio for an experiment may be 
difficult, and this model is therefore impractical. 
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C . Crossf low/T-S Interaction — Finite Amplitude 

For finite-amplitude disturbances, we use a Herbert-type approach to 
model the interaction. Numerical solutions of Eqs. (23)-(26) provide growth 
rate (Eq. (29)), wavenumber, and group velocity-ratio distributions as a 
function of crossflow vortex amplitude. For a chord location of 1.5% and a 
frequency of 325 Hz, these are shown in Figs. 14-16, respectively. The pres- 
ence of crossflow vortices, even with small amplitude, apparently causes 
magnified growth of the oblique T-S waves. The spanwise nonuniformity modu- 
lates the growth rates. In addition, a broad band of wavelengths can be 
observed as a function of amplitude. This is consistent with experimental 
observations for other three-dimensional interactions — specifically, the 
occurrence of three-dimensionality in boundary-layer flows. Figure 16 shows 
the variation of group velocity ratio with amplitude. 

For a chord location of 2.5%, Figures 17 and 18 show growth rate and 
dimensionless wavenumber as functions of frequency and crossflow vortex 
amplitude. Higher-frequency disturbances appear to be more affected by the 
presence of the spanwise nonuniformity, as evidenced by the magnified growth 
rates in Fig. 17. The wavenumber behaves qualitatively the same for the two 
frequencies . 

Figures 19 and 20 show growth rates and amplification factors for a T-S 
wave of frequency 325 Hz and initial amplitude of 0.6%. Results show 
parallel, zero-amplitude and finite-amplitude rates. The growth rates are 
integrated numerically in marching along the trajectory of the wave to find 
amplitude ratio. Figure 20 compares a parallel, zero-amplitude with a finite- 
amplitude model. The results clearly show double-exponential growth for the 
finite-amplitude analysis as compared with the zero-amplitude approach. It is 
clear that the usual transition prediction methods will fail in such a flow; 
premature transition will occur. 

2 .3 Remarks 

The analysis will not be restricted just to the effects of crossflow on 
the T-S instability. The role of T-S waves in affecting the growth of cross- 
flow instabilities will also be examined. Under certain conditions, crossflow 
vortices can also be amplified in the midchord region of hybrid laminar flow 
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airfoils. In this case, we assume the crossflow as the secondary instability 
and proceed with an approach analogous to that above. 

At this time, there is no suitable criterion for establishing transition 
when both crossflow and T-S waves are present. For the lack of anything 
better, a linear relationship is assumed between the n-f actors for crossflow 
transition (n CF ) and T-S transition (n T g)« Usually one assumes 

rtjig = 12 - (1.2) n^ F 

with some error bands, as a transition criterion. The present research will 
be directed toward better understanding of the relationship between n^g 
and . 

The effectiveness of suction will be examined. As was shown by Lekoudis 
(1978) with wall suction, if the percentage of permeable area is small, if 
most of the flow there is directed normal to the surface, if the perforations 
are closely spaced, and if the surface has small permeability, then it is 
acceptable to specify zero boundary conditions for the disturbances. In our 
case, this means that both the crossflow vortices and the T-S waves satisfy 
no-slip and no-penetration. Suction will therefore be introduced into the 
undisturbed mean flow from the Kaups Cebeci code. The stability results of 
this modification of the undisturbed mean flow will be applied to typical LFC 
situations. It is expected that a new capability for analyzing the stability 
of compressible swept-wing flows will result. Whenever appropriate, the work 
will be applied to the presently ongoing transonic LFC experiments being con- 
ducted at NASA-Langley Research Center in the eight-foot Transonic Pressure 
Tunnel (Harvey and Bobbitt, 1982). 
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3. 


Personnel 


The principal investigator for this project is Helen L. Reed, Assistant 
Professor of Mechanical Engineering, Stanford University. William S. Saric, 
Professor of Engineering Science and Mechanics, VPI & SU, has been available 
to consult on this program. The detailed resume of Dr. Reed is contained in 
Appendix C. 

Helen Reed received her Ph.D. in Engineering Mechanics in 1981 from VPI & 
SU and joined the faculty at Stanford University in September, 1982. Before 
coming to VPI & SU, she worked with W. Pfenninger at NASA-Langley on problems 
of supercritical airfoils and laminar flow control. In addition to her inno- 
vative modeling of the suction-strip problem in her dissertation (Reed and 
Nayfeh, 1983), she has developed the state-of-the-art three-dimensional, com- 
pressible stability code (Reed and Nayfeh, 1982) which will be used in the 
proposed work. 

William Saric received his Ph.D. in Mechanics in 1968 from the Illinois 
Institute of Technology. He worked at Sandia Laboratories for ten years in 
the Reentry Vehicle Division and the Atomic and Fluid Physics Division, and he 
has been at VPI & SU for the past eight years. For the past seven years, he 
has been conducting theoretical and experimental studies of boundary-layer 
stability problems for LFC applications (Nayfeh and Saric, 1981; Reynolds and 
Saric, 1982). He is principal investigator on an experimental program of 3-D 
wave interaction in boundary layers that is directly related to the proposed 
work . 

The author is a member of the Heat Transfer and Turbulence Mechanics 
(HTTM) group in the Department of Mechanical Engineering at Stanford, which 
consists of eight faculty members, one research associate, and approximately 
35 Ph.D. candidates. The other members of the HTTM group will also play 
important roles as consultants. 


21 






Fig. 1. X-21 wing coordinate system. 
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Fig. 3. X-21 upper surface suction distribution. 
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Appendix A 
ZEROTH-ORDER PROBLEM 
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NOTE: a: = 0 for crossflow. 

a v , are the chordwise and spanwise components of the wavenumber. 
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